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An infinite array of parallel current-carrying wires is known, from the field of neutral particle 
optics, to produce an exponentially localized magnetic field when the current direction is antiparallel 
in adjacent wires. We show that a finite array of several tens of superconducting Nb nanowires can 
produce a peak magnetic field of 10 mT that decays by a factor of 10 4 over a length scale of 500 nm. 
Such an array is readily manufacturable with current technology, and is compatible with both 
semiconductor and superconducting quantum computer architectures. A series of such arrays can 
be used to individually address single single-spin or flux qubits spaced as little as 100 nm apart, and 
can lead to quantum logic gate times of 5 ns. 

PACS numbers: 



I. INTRODUCTION 



Among the many and varied proposals for constructing 
quantum computers, spintronic solid state devices occupy 
a special place because of the prospects of integration 
with the existing semiconductor technological infrastruc- 
ture 0. At the same time, superconducting devices have 
taken an early lead in demonstrating the viability of the 
building blocks of quantum computing (QC) in the solid 
state, with recent reports of controlled single-qubit oper- 
ations and entanglement generation 

In both the spintronics QC proposals, such as quantum 
dots |3|, U, la, |6( , donor spins in Si p| , electrons on helium 
@ , and the superconducting QC proposals , the ability 
to apply highly localized and inhomogeneous magnetic 
fields would be a definite advantage, if it could be done 
without excessive technical difficulties. In fact the early 
proposals suggested manipulating individual spin qubits 
using such localized magnetic fields, e.g., by a scanning- 
probe tip or by coupling to an auxiliary ferromagnetic 
dot H|, but there are significant speed, controllability 
and other difficulties associated with such methods. Be- 
cause of these difficulties, in particular the spintronics 
requirement to resolve single spins, many alternatives to 
the use of localized magnetic fields have been proposed 
in spin-based QC. These alternatives typically avoid the 
use of magnetic fields altogether: e.g., g-factor engineer- 
ing combined with all electrical control [Tcj . optical spin 
manipulation , or encoding into the states of several 
spins 0, . Other alternatives include gate telepor- 
tation, which requires controllable exchange interactions 
and certain two-spin measurements |14|. and qubits en- 
coded into antiferromagnetic spin clusters, in which case 
the magnetic field needs to be controlled only over the 
length scale of the cluster diameter . In the context of 
superconducting qubits it is also possible to avoid using 
localized magnetic fields by introducing an appropriate 



encoding [l6|. 

Here we revisit the possibility of applying highly lo- 
calized magnetic fields. We show that a scheme inspired 
by magnetic mirrors for cold neutrons [l7| , and more re- 
cently cold atoms |Iglll9ll2Cll2lLl2^ is capable of generat- 
ing a magnetic field that decays exponentially fast over a 
length scale comparable to the spacing between nanofab- 
ricated quantum dots, and has strength and switching 
times that are compatible with QC given available es- 
timates of decoherence times. Our scheme uses arrays 
of parallel current-carrying wires, that is readily imple- 
mentable with currently available nanotechnology Q, |2^| , 
and appears well suited for integration with quantum dot 
nanofabrication methods, as well as with superconduct- 
ing flux qubits and spin-cluster qubits, where the length 
scales are larger. Thus we believe that QC with localized 
magnetic fields deserves a fresh look. 




FIG. 1: Schematic depiction of quantum dots (disks) with 
array of current carrying wires. Array period is a, current 
is / and alternates direction as indicated by arrows. The 
resulting magnetic field profile is drawn schematically. 
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II. EXPONENTIALLY LOCALIZED MAGNETIC 
FIELD FROM AN INFINITE WIRE ARRAY 

In order to have a concrete application in mind we shall 
from now on refer to semiconductor quantum dot spin- 
qubits USUI!. 

However, our results are equally appli- 
cable to other qubits that are manipulated by localized 
magnetic fields, such as superconducting flux qubits [9|- 
The first requirement for single-spin magnetic resolution 
is a magnetic field profile that decays exponentially fast 
over length scales comparable with the inter-spin spac- 



ing. We will now show, in close analogy to results from 
magnetic mirrors, how to produce such an exponentially 
localized magnetic field. The basic design is one of an 
array of parallel current carrying wires, with the current 
direction alternating from wire to wire: see Fig. ^ 

We first consider the idealized case of infinitely long 
wires. In this case the field can be calculated analytically 
(see also ^3,^3)- Let B£ be the magnetic field in the 
a — x, z direction generated by ./V infinitely long wires. 
We add magnetic field contributions from each wire, to 
get the field components from N wire pairs: 
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where I is the current through each wire and k = 2ir/a 
is the reciprocal array constant. This sum can be com- 
puted analytically in the limit N — > oo using the residue 
theorem result 



£ (-!)«/(«) = -7r£^— Res(/,C* 

where Cfe are the poles of /(C), yielding 
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FIG. 2: The magnetic field magnitude from a infinite array 
of infinitely long wires ("doubly infinite"), on a logarithmic 
scale. Note that the field is flat between the wires, and decays 
exponentially in z. 



The z-component result shows the basic point: an expo- 
nentially localized magnetic field can be generated using 
a wire array. The flat top of the sech profile is a useful 
design feature, since it implies no exponential sensitivity 
in the range z < a [f| . The field magnitude 



B°° = |B°°| = V2 [cos(2kx) + cosh(2/b 



-1/2 



(•5) 



oscillates with period a in the x-direction: see Fig. |2] 



III. MULTIPLE ARRAYS 

For the purposes of QC, we should ideally be able to 
address each spin separately. To this end we propose 
to center a separate wire array on each quantum dot. 
Then, as long as the dots are spaced on the order of 
the lattice constant a, we have exponentially sensitive 
addressability of each dot. The introduction of multiple 
arrays is useful in another respect: we can adjust the 
magnitudes and directions of currents in different arrays 
so as to exactly cancel the field at all (or only some) 
other dots except the desired one (or ones). To see this, 
let b(z) = Bf{z)/I = ^-sechkz. The field at position z 
from K arrays of wires, with the jth array having current 
Ij and intersecting the z-axis at position Zj (typically the 
center of one of the dots), is: 



(K-l)/2 
j = -(K-l)/2 



(6) 



Suppose we wish the field to have magnitude Cj at posi- 
tion Zj. Formally, we need to solve: 

B W fe) = Oi, je[-(K-l)/2,(K-l)/2}. (7) 
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magnetic field (units of ;i //a) 




FIG. 3: Magnetic field generated by a single doubly infinite 
wire array centered at z — (solid), vs the field generated by 
K = 5 such arrays (dashed), with currents chosen to cancel 
the field at positions z = ±2o, ±o. 

This is a linear system of K equations in the K unknowns 
Ij , so it can always be solved in terms of the K positions 
Zj. E.g., the fields with and without the correction are 
shown, for K = 5, in Fig. [S] 

IV. FINITE SIZE EFFECTS 

For a finite system (TV < oo wires, finite length and 
thickness, non- ideal shape, etc.), we can only expect the 
above results to hold to an approximation. Much of the 
theory of corrections to finite size effects has already been 
worked out in |l9l l22j |. in the context of atomic mirrors. 
Since for atomic mirrors the primary concern is specular 
reflection, there the focus was on reducing the variation 
of the field magnitude in the planes parallel to the wire 
arrays. For us this criterion is unimportant; instead, our 
focus is on making the field as localized as possible along 
the z-axis. 

The first important conclusion in the case of a finite 
number of wires N is that there is a transition from 
exponential to quadratic (i. e., 1/z 2 ) decay. [2^ The 
transition takes place at the inflection point (zero second 
derivative) of B N (z); however, this is difficult to obtain 
analytically To roughly estimate the transition point we 
compute where the field from a single pair of wires, posi- 
tioned at the edge of an array of N wires (y — ±7Va/4), 
generates a field of magnitude equal to that from an in- 
finite array: 

Hal Na/4 = /x J 27rz/a , . 

7T z 2 + (Na/A) 2 a ' 1 ' 

Since the transition happens for x / a <C N we neglect x/a 
in the denominator. The solution is then 

ira , ,ttN . ,„. 

z t K, — ln( ) (9) 

20 V 4 ' y ' 

Numerical calculations show that Eq. overestimates 
the position of the transition point by about a factor of 
5; however, after this correction is made, analytics and 
numerics agree well across several orders of magnitude of 
N. 



This logarithmic dependence might appear to pose a 
severe scalability constraint on our method. However, 
this is not the case when we take into account the thresh- 
old for fault tolerant quantum error correction |24j]. For, 
it follows from the threshold result that we only need to 
make the ratio of the peak field (applied to the desired 
spin) to the residual field smaller than, say I0 -4 . The 
crucial question thus becomes for what N this can be 
achieved, and this brings us to the idea of "endcaps" . 

As observed in [l9ll2lll22]| . near the center of the array 
the magnetic field that would be produced by the semi- 
infinite array of "missing" wires is the same (to first order 
in 4/iV) as that of a pair of wires carrying current J cap = 
J/2 and placed with their centers shifted by a/4 from 
the outermost wires in the array. Thus, to cancel this 
field, one can simply place two "endcap" wires carrying 
currents ±1/2 at these positions. In the context of atomic 
mirrors this is important to improve flatness, and hence 
specular reflection. In |23 it was observed that flatness 
can be further enhanced by using an odd number of wires. 

Some of these schemes can be used to improve field 
localization, a criterion not considered orginally. For in- 
stance, we find that the number of wires can be reduced 
drastically - from TV w 1 4 to as few as N = 22 wires - 
when endcaps are used to achieve a residual field smaller 
than I0 -4 . By contrast, using an odd N is disastrous for 
localization: for example, residual fields appear at the 
3% level for N = 23, even including endcap correction. 
Higher (even) wire number increases the robustness of 
the cancellation against experimental uncertainty in the 
current and position of the endcaps: N = 30 is required 
to maintain B/hqI < I0~ 4 for a fractional current vari- 
ation of ±10 -3 and a positional uncertainty of ±2.5 nm, 
as shown in Fig. 0] Finally, we note that the corrections 
arising from the finite length of the wires and the short, 
perpendicular connecting wires, can also be compensated 
for by the use of judiciously placed compensating wires 



V. FEASIBILITY AND IMPLEMENTATION 
CONSIDERATIONS 

We now come to estimates of whether the fields and 
size scales required are feasible in practice. Let us first 
calculate the magnetic field strength required for single- 
qubit operations. A spin can be rotated by a relative 
angle (j> = g^sSr /2h by turning on the field B for a 
time t (where g « I is the (/-factor, and fj,B is the 
Bohr magneton). Recent estimates of dephasing times 
are 50 fis for electron spins in GaAs quantum dots (a cal- 
culation, assuming spectral diffusion is dominant) |25j . 
and a measurement of 60 ms for T2 of phosphorus donors 
in Si j2|| . If we use the more pessimistic of these num- 
bers jand assume a fault tolerance threshold of I0~ 4 for 
QC [2J|, we may estimate the desired operation time as 
t ~ I0~ 4 x 50 /is — 5 ns for an angle (f> = ir/2. Thus the 
desired field strength is B = 2h(f>/gfiBT ~ 7mT. 
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depth, the critical current density is due to depairing: 
j c = (2/3) 3 ^ 2 B c /(p n X) 29J . These constraints are com- 
patible when A < £o, i.e., mostly type I superconductors. 
For Nb, B cl = 0.206 T, A = 52 nm, and fn = 39nm |30|. 
so j c = 1.7 x 10 12 A/m 2 . Note that high-temperature 
superconductors typically have higher A and thus lower 
j c . In any case, the critical current density of Nb is more 
than is required: a gate time of r = 5 ns would require 
j » 0.1 Jc . 



FIG. 4: Deviation |B°° - B x \ of the field of a finite array 
(solid lines) from the field |B°°| of an infinite array (grey 
dashed line, shown for comparison). This array has N = 30 
wires (including two endcap wires), such that main wires are 
at ±a/4, ±3a/4, . . . ± 27a/4, and endcap wires are at ±7a, 
carrying current 7 cap = 7/2. Field deviations are plotted 
for a range of endcap currents (±10 _3 7 cap ) and positions 
(±a/100 = 2.5 nm). Note that the residual field never ex- 
ceeds the threshold of 10 -4 . 

To evaluate the feasibility of such a specification, we 
consider an array with periodicity a = 250 nm, wire ra- 
dius r = 50 nm, N = 32 wires, and length L y — 10 /im 
along the y-direction. These dimensions are compati- 
ble with the 100 nm length scales of quantum dots [6j. 
In order to reach the desired field strength of 7mT, 
I = aB^° j hq w 1.4 mA would be required. However, 
decoherence due to heating with such a current could 
be a major issue. An upper bound estimate for the 
required temperature T can be given by T <C Ez/ks, 
where Ez = gfisB/2 is the Zeeman splitting of the 
spins in the applied magnetic field B. In our case, we 
have constrained B by the gate time r, so we can write 
k B T < nh/2T, or T < 2.4 mK for r = 5 ns. This is 
feasible with dilution refrigeration technology if the heat 
load is on the order of 100 pW [27]], comparable to the 
dissipation of quantum dots 0] . 

For normal metal wires, such a heat load restriction 
is prohibitive. The power P dissipated is P = j 2 pA£, 
where j is the current density, t = N(2L y + a) is the 
total wire length, and p is the resistivity. Below 10 K, 
oxygen-free copper can have p«3x 10~ n fim. |28( At 
P < 100 pW, I < 3.6 juA, which would give r > 2.0 ^s 
- nearly a thousand times slower than our original goal. 
Although this would be acce pta bly fast if the decoher- 
ence time were 60 ms, as in |26J, the high-purity p we 
have used is optimistic for nanofabricated wires, and the 
heat generated would increase linearly with the number 
of qubits manipulated. 

One can circumvent resistive heating by using super- 
conductors. A wire with radius r < £bj where £o is the 
coherence length, can also avoid heating mechanisms as- 
sociated with vortex movement through the supercon- 
ductor. For r~Aor smaller, where A is the penetration 



VI. CONCLUSIONS AND EXTRAPOLATIONS 

In conclusion, our results indicate that QC with local- 
ized magnetic fields deserves renewed consideration. We 
have shown that a method to produce exponentially de- 
caying magnetic fields using an array of current-carrying 
wires, known in the cold neutron and atom optics com- 
munities, is adaptable to solid state quantum computer 
implementations. Our estimates indicate that in all re- 
spects the method is technologically feasible, provided 
superconducting wires with sufficiently high critical cur- 
rent density (such as Nb) are used. 

Our work is motivated by the quest to perform single- 
spin or flux-qubit rotations, which is a component of 
a universal set of quantum logic gates. The geometry 
shown in Fig. ^ yields a field that is localized in the z- 
direction; in order to perform arbitrary single-qubit ro- 
tations we need to localize the field along another, per- 
pendicular direction. An independent field vector can 
be produced by a second set of interleaved arrays placed 
at 45° with respect to the original arrays, with current 
flowing along (z + y)/ \/2. With a judicious array place- 
ment and na = d/ v2, for qubit spacing d and any integer 
n, the field direction will be along x at all qubits (this 
design will need to be optimized similarly to our consid- 
erations above - an issue we do not intend to address 
here). The additional spatial constraints would require 
only a four-fold increase in a for the same currents and 
wire sizes. If on the other hand, introducing a second ar- 
ray is undesirable, "software" solutions using recoupling 
and encoding techniques have been developed to still al- 
low for universal QC These techniques would be 
considerably simplified by the ability to perform single- 
qubit operations along one direction. 
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